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Abstract 

We consider systems which exhibit typical critical dependence of the 
specific heat: Ac a (Tc - T)"^ (T < Tc); Ac a (T - Tc)"^' (T > Tc) 
where 7,7' are critical exponents (7 = for Ac = Acp,jv, 7 = for 
Ac = Acv.iv), as well as, the case when Ac a (In | Tc — T |)" (a = i, 
uniaxial ferroelectrics; a — 1, liquid He^). Starting from the critical 
behaviour of the specific heat we find the Gibbs (Helmholtz) potential in 
the vicinity of the critical point for each case separately. We derive in 
this way many exact critical relations in the limit T —> Tc which remain 
the same for each considered case. They define a new class of universal 
critical relations independent from the underlying microscopic mechanism 
and the symmetry of these systems. The derived relations are valid for a 
very broad class of magnetic, ferroelectric and superconducting materials, 
as well as, for liquid _ffe*. 

Subject classification: 05.70.Fh, 05.70.Jk, 74.25.Bt, 67 



1 Introduction 

Critical phenomena in solids, connected with second order phase transitions, 
belong to very attractive and intriguing problems in physics. Many theoretical 
models and methods have been applied to resolve these interesting questions 
(cf e.g. Refs [l]-[5] for a review and the original papers cited therein), how- 
ever, the most important, general and successfuU method, applied in this area 
of physics, seems to be still the phenomenological approach, initiated by the 
pioneer works of Landau, Refs [6], [7]. This approach has succsessfully been 
applied to superconducting systems resulting in many useful critical relations, 
derived in Ref. [8]. For many other systems (magnetic, ferroelectric, liquid 
He^, etc.) this approach unfortunately fails because of critical exponents (cf 
e.g. Refs [l]-[5]). We can, however, improve this approach and calculate the 
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Gibbs (Helmholtz) potential near the critical point from the observed specific 
heat anomaly [T <Tc;T >Tc). Similarly to Ref. [8] we include into the con- 
siderations the critical behaviour of the chemical potential as indicator of phase 
transitions in a more general sense (see Ref. [8] and papers cited therein). We 
derive in this way many interesting thermodynamical relations (cross-relations) 
at the critical point in the limit T — > Tp. These exact relations seem to be of 
fundamental importance for critical phenomena in solids and should generally 
be valid for a very broad class of magnetic or ferroelectric systems, as well as, 
for superconducting systems as a special case. Surprisingly, they are also valid 
in the case of logarithmic divergence of the specific heat (uniaxial ferroelectrics 
and liquid He^). The derived relations possess very similar form to the critical 
relations presented in Ref. [8] and define a new class of exact relations at the 
critical point independent from the underlying microscopic mechanism leading 
to second order phase transitions and independent from the symmetry of the 
system. 

2 Phenomenological Gibbs (Helmholtz) poten- 
tials and critical relations 

Performing experimental measurements of the temperature dependence of the 
specific heat Cp^N at constant pressure p and constant number of particles N 
we can find that for magnetic, ferroelectric or superconducting systems Cp,jv ex- 
hibits a pronounced anomaly at the critical point T = Tc (cf e.g. Refs [l]-[5]). 
Independently from the considered system (here: magnetic, ferroelectric or su- 
perconducting) it is convenient to introduce a reference system which possesses 
exactly the same structure (symmetry), atomic masses, etc., the specific heat 
of this system c^^]^ exhibits, however, no anomaly at T = Tc- Let us introduce 
an auxiliary quantity 

AcpMT)^CpMT)~cl%{T). (1) 

A very broad class of materials (magnetic, ferroelectric or superconducting) 
shows a general, critical behaviour (see e.g. [l]-[5]) of the type^ 

where a, a' are critical exponents, depending on the material (a, a' > 0). For 
magnetic and ferroelectric systems Acp,jv usually diverges at the critical point, 
for superconducting ones Acp,iv exhibits only a finite jump. For simplicity, we 

^The case of the logarithmic divergence of the specific heat will be considered in the next 
Section. 
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restrict ourselves to the case T < Tc and we assume (see expression (2)) that 
in the vicinity of the critical point 



Acp,N(T) - A{p,N)(Tc - T)-« (3) 

where A{p, N) is a parameter. The critical temperature Tc should also be con- 
sidered as a function ofp and N {Tc = Tc{p, N)), as suggested by experimental 
results. The thermodynamics of the system in the vicinity of the critical point 
(r < Tc) can be described by the Gibbs potential Z{T,p,N). We can formally 
write 

Z{T,p,N) = Zo{T,p,N) + AZ (4) 

where 

AZ = Z- Zo, (5) 

Zq is the Gibbs potential of the reference system and c^"]^ = — T ^ g^j' ^ 

The analytical form of AZ should necessarily be chosen in such a way that the 
thermodynamical relation 

reproduces the critical behaviour of Acp^N, given by the formula (3). It is easy 
to find that the leading term of AZ should have the general form^ 

AZ = D{p, N){Tc - T)^-" (7) 

where D{p, N) is a parameter. Really, when applying (6) and (7) we obtain 

Acp,jv = -T{2-a){l-a)Dip,N)iTc-T)-" 
« -Tc{2-a)(l^a)D{p,N){Tc-T)-'' 
= A{p,N){Tc-T)-^ (8) 

where 

A{p,N) = -Tc{2-a){l-a)D{p,N). (9) 

Thus, the asymptotic form of the Gibbs potential (4) in the vicinity of the 
critical point [T <Tc) has to be equal to 

Z{T,p, N) = Zo{T,p, N) + D{p, N){Tc - Tf-^. (10) 



^It is possibly to find exactly when integrate the expression (6) with the use of the 
formula (3). The result is given in the Appendix A. To avoid the presentation of very long 
expressions we restrict ourselves to consider only the leading term in AZ, given by the formula 
(7). This approach (see Appendix A) leads to the same, exact critical relations (18)-(22) as 
in the case when the complete expression for AZ (see the formula (A. 6)) is used. 
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The phenomenological relations (3) and (10), resulting from the experimental, 
critical behaviour of a broad class of systems can easily be used to derive many 
critical relations. The chemical potential /i = (|^)t.p (pq = (^ff)T,p) and the 
volume of the system V = (^)t,jv {Vq = {^§^)t,n) can easily be found from 
(10). We obtain 



p 



and 



(2-a)D(Tc-T)i-"(^^) (11) 



N 



fdTi 



+ (2-a)D(Tc-T)i-" — ^ (12) 



V dp 



N 



where we have taken into account that Tc = Tc{p,N). With the use of (11) 
and (12) we obtain 



-(Tc-T)-«(2-a)[(|^) (Tc-T) 



(13) 



V dN 



p 



^ = -(^c-T)-"(2-a)[(^ ) (Tc-T) 



N 



(14) 



+(2-a)(Tc-r) 



dD\ fdTc 



ONJ^ V dp 



)] 
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{dNj 



(Tc-T)- 



+ (2-a)(Tc 



T) 



+{2-a){l-a)D 



and 



) = (rc-T)-"{( 

/ T,N ^ \ 



+(2 - a){Tc - T) 
+{2-a){l-a)D 



dD 
dN 

dTc 
dN 



{Tc 

N 

dD 
dp 
dTc 
dp 



N 



Tf 

dTc 
dN 



Tf 

dTc 
dp 



D- 



d^Tc 

dN"^ 



(16) 



+ D 



N 



d^Tc 
dp^ 



N 



N 



(17) 



where AX = X{T) — Xo{T) has exactly the same interpretation as in the 
formula (1). It is interesting to note that for magnetic (ferroelectric) svstcnns 
Acp.jv, A(|^)p,jv, ^i§^)p,N, A(|^)t,jv, A(|^)t,p and M^W^n , generally, 
diverge at the critical point. However, when e.g. divide (13)-(17) by (3) the 
"dangerous" term {Tc — T)~" cancels out (see also (9)) and taking the limit 
T — » T^ we obtain the following relations 

Mw)p,N _ fdlnTc' 



lim 



Ac 



'p.N 



lim 



lim 



lim 

Adf) 



A(f^)p,jv 



Ac 



'p,N 



T,N 



ACp^N 



and 



lim 



T.N 



Ac. 



p.N 



dN 

dlnTc 
dp 

dlnTc 
dN 

dlnTc 
dN 

dlnTc 
dp 



N 



[dTc 



dp 

dTc 
dN 



N 



N 



dTc 



c 



dp 



(18) 



(19) 



(20) 



(21) 



(22) 



N 



It, however, means that the limits exist and the derived relations are exact'^. 
For classic superconductors {a = 0) all the quantities (3) and (13)-(17) are finite 

■^Thc additional terms to AZ (sec Appendix A) modify only the expressions (11)-(17). By 
forming the quotients (18)-(22) in the limit T their contribution is equal to zero (see 

Appendix A) and therefore the critical relations (18)-(22) are, in fact, exact. 
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and therefore we can drop the lim in the expressions (18)- (22) and we obtain 

the formulae (29)-(33) from Ref. [8]. In orther words the formulae (18)-(22) are 
valid for a broad class of magnetic, ferroelectric and superconducting systems. 



Because of Maxwell relations (f = - (|^)^_^, {§§)^^^ = - 
and (f^)j.p — (if) there are only five independent relations of the type 
(13)-(17). Because A (f = -A (fl)^^, A (|f )^^ = -A (|^)^^^^ and 
A {§^)rpj^ ~ A (ip) corresponding critical relation with the presence 

of A (if )^^, A (w)t,p and A (|^)^pCan easily be obtained from (19), (18) 

and (20), respectively. Besides, it is very easy to obtain many other relations 
when e.g. divide (14) by (13) and take the limit 

lin, A(^W^Ui^ (23) 

and so on. Therefore there is no need to write all of them explicitely here. It 
is also possible to calculate the derivatives (^^)p , i^§^)N from the relations 
(18), (19) and isert them into (20)-(22). In such a way we obtain cross-relations 

lim 7 = -Tc[ lim • [ T^""] (24) 



lim =-Tc[ lim ^pZP^]^ (25) 



and 



hm A!^ = _Te[ lim i^^^-^ll (26) 

In the case of superconducting systems we can drop the lim and the formulae 

(24)-(26) coincide with the formulae (34)-(36) from Ref. [8]. There are, however, 
again many other possibilities to form similar cross-relations which can be found 
with ease when starting from the formulae (3), (9) and (13)-(17). 

We can also derive additional relations for systems at constant volume V. 
Similarly to (2) we can write (cf e.g. [l]-[5]) 

Ac.,.(T)a{ nil (27) 
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where a and a' are critical exponents. We restrict ourselves again to the vicinity 
of the critical point for T <Tc and we can assume that 

Acv,N = A{V, N)iTc - T)-« (28) 

where A(V, N) is a parameter. The corresponding leading term of the Helmholtz 
free energy has to be in the form (cf also (10)) 

F{T, V, N) = Fo{T, V, N) + D{V, N){Tc - Tf'" (29) 

with D{V,N) as a parameter and Tc = Tc{V,N). Here Fq can be interpreted 
as the Helmholtz potential of the reference system. Repeating the calculations 
with the use of (29), quite similar to the presented above with the use of the 
Gibbs potential (cf (11)-(17)), we can find the exact relations similar to (18)- 
(22)^. We obtain 

M^)v,N _ fdlnTc 



lim -oi^','- ^ (30) 
^{^)v,N _ fdlnTc 



T^T- Acv,N \ dV Jj^' 

lim = - ('-^) C-^) , (32) 



Acv,N V dN Jy\dV 

^^-i'-^)m, (33) 



and 



Acy,iv V dN J,,\dN 



^{w)t,n fdlnTc\ (dTc 



There are again only five exact expressions of this type because of the Maxwell 
relations (M)^_^ = - (t)^^, (||)^,^ = " and (^)^ 



V,N / V,N V^'/T.V 

. There exists also a freedom to form another quotiens what re- 

T.N 



suits in many other expressions easy to obtain (we don't write them explicitely 
here). The analogous cross-relations, similar to (24)- (26) have the form 



— = Tc hm '— • hm — '- 

A.Cv,N T^T- ACv,N T^T- ACy^N 



^Thc same critical relations (30)-(34) can also be obtained with the use of the exact 
Helmholtz potential, similarly to systems at constant pressure p (see Appendix A) and there- 
fore the formulae (30)-(34) are exact. 
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\dN J 



= -MJ^. ^T^]^ (36) 



and 



lim ^P^=Tc[ lim 4^]^ (37) 



The other cross-relations arc also very easy to find. In the case of supercon- 
ducting systems we can drop the lim in these expressions and we obtain the 

T^T- 

formulae (46)- (53) from Ref. [8]. Similar calculations can also be repeated in 
the vicinity of the critical point for the case T > Tc, using the form of the spe- 
cific heat for this case (see (2) and (27)). It leads, however, to the same critical 
expressions (18)-(22), (24)-(26) and (30)-(37) where the Hmit T ^ should 
be replaced by T — > . In other words, we can replace in all these expressions 
the limit T —> by a more general form lim . 

T—^Tc 



3 Logarithmic anomaly of the specific heat 

It is interesting to see whether the derived critical relations (18)-(22), (24)-(26) 
and (30)- (37) are also valid in the case of logarithmic divergence of the specific 
heat, when 

Acp,;v(T)cx(ln(Tc-T))« (38) 

where a | for uniaxial ferroelectrics (see e.g. Ref. [5], p. 384) or a = 1 for 
liquid He* (see e.g. Refs [1], [9])^ 

For simplicity, we restrict ourselves to the case T <Tc where 

Acp,A,(T) = Aip, N)[ln{Tc - T)^ (39) 

and A{p, N) is a parameter. It is easy to see that the leading term of the Gibbs 
potential, deduced from the expressions (6) and (39) should have the form^ 

Z{T,p,N) = Zo{T,p,N) - ^D{p,N){Tc - TfMTc - T)]« (40) 

^According to different authors for liquid He'^ the formula (2) can also be valid for a, a' 
very small (cf e.g. [3] and papers cited therein). Thus, the critial relations, derived in the 
preceding Section are autliomatically valid in this case. 

''The exact form for AZ is given in Appendix B. Similarly to the preceding Section, the 
additional terms in AZ does not contribute to the final relations (18)-(22) which remain 
exactly the same. Because of extremally long expressions we present here only the calculations 
resulting from the leading term in the formula (40). 
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where D[p, N) is a parameter. Applying the formulae (5) and (6) we obtain 
Acp,N{T) = TD{p,N){[\n{Tc-T)r 

+ ^[ln(Tc - T)]«-i + ^^^[ln{Tc - T)]-^}. (41) 

In the vicinity of the critical point the second and third term in the parenthesis 
of (41) can completely be neglected (a = ^, 1) in comparison to the first, singular 
term. Thus, wc can write 

Acj,,iv(r) « TD{p, N)[\n{Tc - T)]" « TcD{p, N)[\n{Tc - T)]«. (42) 

In other words the formula for the Gibbs potential (40) reproduces the critical 
behaviour of the specific heat (39) when we assume that 

A{p,N)=TcD{p,N). (43) 

The chemical potential = (f^)T,p (mo = (^^)t,p) and the volume of the 
system V = {^)t,n {Vq = {^)t,n) can be found from (40). We obtain 

= ^'-\{^) iTc-Tf[HTc-T)r 



and 



D{Tc - T){MTc - TT + 2 [ln(Tc - T)^-'} [-^ j (44) 



V = Vo-\[^^ ^{Tc-TnHTc-T)]- 
-D{Tc - T){[\n{Tc - T)r + ^[ln(Tc - T)]«-i} ] (45) 

N 



2' ' ^ " ' \ dp 



where we have taken into account that Tc = Tc{p,N). It is easy to see that 
according to (44) and (45) we find 



(dTc\ 



D{[ln(Tc-T)]« 



p 



+^[HTc - T)]«-i + ^i^[ln(rc - T)]«-^}, (46) 
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'dir 
- dp . 



N 



{Tc - T){[\n{Tc - T)Y + -[ln(Tc - T)]«-i} 



+ 



(t) 



D{[\n{Tc - T)f 



N 



+ ^[ln(Tc - T)f-' + ^^^[ln(Tc - T)]-^}, (47) 



-2i^::^r)iTc-TnHTc-T)r 



\dpdN) 



dTc 
dp 



{Tc-T){[ln{Tc-T)]'' 



N 



+ -[ln{Tc-T)r-'} 

+ |[ln(Tc - T)r-' + ^^^[ln(Tc - T)]"-^}, (48) 



1 fd'^D\ 
+ ^[ln(rc-T)]-i} 

-^[( W) ^'^tM^C - T)]" + |[ln(Tc - T)]«-i 

+ ^^^[ln(Tc-T)]-2| (49) 
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and 




) 



T,N 




+ ±^[lr,(Tc-T)r-'} 



a-l 



(50) 



where, as before, AX = X{T)~ Xo{T). It is easy to see that dividing (46)- (50) 
by (39) with the use of (43) and taking the limit T we obtain exactly 

the same critical relations (18)-(22) and also (24)-(26). The calculations in the 
case T > Tc are quite similar to the case T < Tc and we do not repeat them 
here. They also lead to the same results (18)-(22) and (24)-(26) where the limit 
T should be replaced by T ^ Tq . In the other words, similarly to the 

preceding Section, we can simply use the limit ^lim in these expressions. The 

case of the system at constant volume V can be treated parallel to the lines 

dcsribcd above and in the preceding Section. It, however, leads again to exactly 
the same critical relations (30)- (37). Thus, we see that in the case of the specific 
heat anomally (2) and (38) the same critical relations are valid. 

4 Conclusions 

Using the phenomenological approach to second order phase transitions we have 
shown that independent from the type of the considered system the critical re- 
lations (cross-relations) derived in this paper are exactly the same for all of 
them. We have shown in this way that quite different systems possess the same 
universal critical behaviour of the quotients. It, however, means that the exact 
critical relations, derived in this paper, define a new class of universal, critical 
relations valid for a very large number of magnetic, ferroelectric or supercon- 
ducting systems including also uniaxial ferroelectrics and liquid He^. 

5 Appendix A 

Because of the thermodynamic relation 




(A. 1) 
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we can write 



AS = I ^^dT + Ci ip, N) (A. 2) 

where Ci{p, N) is an arbitrary constant. In the vinicity of the critical point the 
condition <C 1 is fulfilled and therefore we can write 

T Tc-{Tc-T) Tc 1 - (^^-^) 7bf-f ■ ^ ' ' 

Tc (=U ^ 

Thus, we obtain (see the formula (3)) 

A5 =^Mg J^J(Tc-Ty-«dT + Ci(p,7V) 



(A. 4) 



Because a <1 the first term in (A. 4) is going to zero at T = Tq- At the critical 
point AS* = Q {S = Sq) and therefore the arbitrary constant Ci(p, A'') has to be 
equal to zero. 

Starting from the expression (A. 4) we obtain 

AZ = - J ASdT + C2ip,N) 

At n.r\ °° IT T"i'-<»+2 

At the critical point, however, the first term in (A. 5) vanishes. Because at the 

critical point AZ = {Z ~ Zq) the arbitrary constant C2{p, N) has to be equal 
to zero. Thus, we can write the exact expression for AZ in the form 

A 7 _ A(p,Af) {Tc-Tf-" 
- (2-a)(l-a) 

(A. 6) 

-A{p,N)Y. 

/=1 [l-a+2){l-a+l) 

and we see that the first term (leading term) in (A. 6) has the form (7) (see also 
(9)). This term properly reproduces the critical behaviour of the specific heat 
(sec formula (8)) and therefore we have neglected in further considerations (see 
(11)-(17)) the second term in (A. 6). It is, of course, possible to take into account 
the full expression (A. 6) for AZ but it is completely irrelevant. The presence 
of the second term in (A. 6) modifies the relations (11)-(17) producing very long 
expresions. It has, however, absolutely no influence on the critical relations 
(18)-(22) because the additional terms in (11)-(17) make no contributions in 
the limit T — > when forming the quotients (18)-(22). 
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6 Appendix B 

Starting from the expression (39) and using (A. 3) we find 



AS =1 i^^dT + Ci {p, N) 

(B. 1) 

= [ln(Tc - T)r E %^Wi{Tc -T) + C,{p, N) 

where 

oo 

Wi{Tc -T)=J2 (^i,m MTc - T)]"" , (B. 2) 

m=0 

ai,o = l, au = -ifT. — a,,„ = (-l)- °<°-^\— 5°-'"+^^ (B. 3) 

and Ci{p, N) is an arbitrary constant. The first term in (B.l) vanishes when 
T ^ (a = Because AS" should vanish at T = Tc (S = Sq) the 

arbitrary constant Ci{p, N) should also be zero. Finaly we find that the exact 
expression for AZ has the form 

AZ =- (Tc - Tf [ln(Tc -T)fx 

oo , ,j oo f oo "1 

X E tS^+^W) I MTc - T)]-^ E bi,m,n MTc - T)]-" I 

+C2{P,N) 



(B. 4) 



where 



I _ 1 L _ a-m I, _ f 1 ^l^ a(a-m)(a-m-l)-....-(a-m-K+l) 

(B. 5) 

and C2{p,N) is an arbitrary constant. Because the first term in (B.4) vanishes 
in the limit T the arbitrary constant C2 (p, N) should be equal to zero 

{AZ = {Z = Zo) at T = Tc). The expression (B.4) can be rewritten in the 
short form to be 



AZ = - ^^P^^\ tc - Tf MTc - T)]" {1 + R{Tc - T)} (B. 6) 

where R{Tc — T) contains the terms with l,m,n= 1, 2, ... (cf (B.3) and (B.5)). 
The first (leading) term in (B.6) has the form of the second term in expression 
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(40) where D{p, N) is related to A{j), N) by the formula (43). It is, in principle, 
possible to perform the calculation (formulae (44)- (50)) starting from the full 
expression (B.6). It, however, leads to extremely long formulae and therefore 
we have presented the results using only the leading term. The terms, generated 
by R{Tc — T) in (B.6) produces additionally extremely long expressions in (44)- 
(50) which arc completely irrelevant when forming the quotients (18)-(22) and 
taking the limit T because their contribution to the expressions (18)-(22) 

is equal to zero in this limit. Therefore the critical relation (18)-(22) are also 
exact in the case of the logarithmic divergence of the specific heat (39). 
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